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ABSTRACT 

Using Chandra observations we have measured the energy-resolved dust-scattered X-ray halo around 
the low-mass X-ray binary GX5-1, which shows signs of both singly- and multiply-scattered X-rays. 
We compared the observed X-ray halo at various energies to predictions from a range of dust models. 
These fits used both smoothly-distributed dust as well as dust in clumped clouds, with CO and 21 cm 
observations helping to de termine the position of the clou ds along the line of sight. We found that 
the BARE-GR-B model of lZubko. Dwek fc Arendt I (|2004Tl generally led to the best results, although 
inadequacies in both the overall model and the d a ta lim it our conclusions. We did find that the 
composite dust models of lZubko. Dwek fc ArendTI l)2004|) . especially the "no carbon" models, gave 
uniformly poor results. Although models using cloud positions and densities derived naively from CO 
and 21 cm data gave generally poor results, plausible adjustments to the distance of the largest cloud 
and the mass of a cloud in the expanding 3 kpc Arm lead to significantly improved fits. We suggest 
that combining X-ray halo, CO, and 21 cm observations will be a fruitful method to improve our 
understanding of both the gas and dust phases of the interstellar medium. 
Subject headings: dust — scattering — X-rays: binaries — X-rays: ISM 



1. INTRODUCTION 

Although interstellar (IS) dust is a significant compo- 
nent of dense molecular clouds, measuring the properties 
of these grains is difficult. Except in the outskirts of the 
clouds, the inherently large extinction prevents optical 
and UV measurements except in unusual circumstances. 
As a result, information about the size distribution and 
composition of grains in dense clouds has by necessity 
been extrapolated from observations of less-dense clouds 
or modeled based on IR measurements. We show that X- 
ray scattering observations, combined with CO spectral 
line measurements, can put useful limits on the allowed 
dust models. 

Our X-ray source is the low-mass X-ray binary 
(LMXB) G X5-1 at (I b) = (5.0 77,-1.019). GX5-1 is a Z 
source (see Ivan der Klis 1 119951 and references therein), 
and is the second brightest persistent Galactic X-ray 
source after the C rab Nebula. The se condary has been 
imaged in the IR ( Jonker^et^aL^ 2000) but the IR spec- 
tra measu red bvlBand vopadh vay et al. 1 l|2008D could not 
be typed. iChristian fc Swank I <|1997ir estimated the dis- 
tance to be 9 kpc, although they warned that this is 
probably an upper limit to the distance and assigned an 
error of 2.7 kpc to the value. Using X-ray spectra taken 
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by the Einstein satellite they also measured the column 
density to be N H = 2.54 ± 0.19 x 10 22 cm" 2 , although 
this re sult depends o n the assumed spectral model. Re- 
cently, lUeda et al. I p005|) used a Chandra HETG ob- 
servation to measure the depth of absorption edges of a 
number of IS metals along the line of sight to GX5-1. 
Based on these measurements, they determined the to- 
tal A^h^ 2.8l fg x 10 22 cm~ 2 ; this result has larger error 
bars than the IChristian fc Swank I l|1997j) value, primar- 
ily because it does not assume a spectral model a priori. 

The X-ray halo around GX5-1 has been observed 
by many X-ray satellites, albeit with lower angular 
or energy resolution than is pos sible with the Chan- 
dra ACIS. iPredehl et al. I l)1992f) used a lunar occul- 
tation of GX5-1 to measure the scattered halo di- 
rectly, finding a total halo intensity (in the ROSAT 
band) 28% of the source intensity . GX5- 1 was also in- 
cluded in the IPredehl fc Schmitt I l|l995f) survey of X- 
ray halos detected with ROSAT, where they measured 
a halo intensity 25.4% of the source intensity. Recently, 
iXiang. Zhang, fc Yao I ^00^1 analyzed 17 X -ray sources, 
including GX5-1, using a new technique based on the 
zero-order image of Chandra HETG data. They found a 
lower total halo intensity (9.9%), as expected since the 
Chandra bandpass extends to higher energies where halo 
intensities are smaller. They also measured a total col- 
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umn density Nr= 2.0 x 10 22 cm -2 and fit their data to 
determine the relative dust positions along the line of 
sight, finding that 90% of the dust was within 10% of 
the source. 

2. THEORETICAL FRAMEWORK 

lOverbeckl (^965) first described astrophysical X-ray 
halos, which result from small-angle scattering of X-rays 
inside solid matter. The technique has long been used 
in laboratory studies of virii or proteins in solution, as 
careful measurements of the radial profile can show the 
densi ty and shape of the particles l|Guinier fc Fournetl 
1955). The theory of X-ray halos in an astrophysical 
settin g has been described in detail by a number of au- 
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the theory here. 

The fundamental quantity is the differential scattering 
cross section da/dfl, which can be calculated using either 
the exact Mie s olution or the Rayleig h-Gans (RG) ap- 
proximation; see ISmith fc Dwe"k1 1)1998[) for a discussion. 
As the RG approximation is analytic it demonstrates how 
the various parameters scale: 
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where a is the grain radius, Z the mean atomic charge, 
M the mean atomic weight (in amu), p the mass den- 
sity, E the X- ray energy in keV, F(E) the atomic scat- 
tering factor ijHenke l fl98l\ (9 sca the sc attering angle, 
and 4> 2 (0 sca ) the scattering form factor (Mathis & Lee 
119911) . For homogeneous spherical grains, the form fac- 
tor is given by <& 2 (0 sca ) = 3(sinu — u cos u)/u 3 where 
u = 47rasin(fl sra /2 ) /A ~ 2ira9 !irs ,E /he. 

ISmith fc Dwek I l|1998j) showed that the RG approxi- 
mation will overestimate the total scattering if the en- 
ergy of the X-rays (in keV) is not substantially larger 
than the size of the dust grains (in pva), and suggested 
2 keV as a minimum energy for most ISM dust models. 

By integrating the scattering cross section over the line 
of sight geometry, the source spectrum, and the dust size 
distribution we get (considering single scatterings only) 
the halo surface brightness at angle 6h from the source: 

I S cM) = FxN H f dE S(E) f da n(a) x (2) 
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where Fx is the total source flux, JVh is the hydrogen 
column density, S(E) is the (normalized) X-ray spec- 
trum, n(a)da is the dust grain size distribution, and 
Osca ~ thetah/{l — x). Here f(x) is the density of hydro- 
gen at distance xD from the observer divided by the line 
of sight average density, wh ere D is the distance to the 
source ijMathis fc Lee 1119911) . 

2.1. Multiple Scattering 

If the column density is sufficiently l arge, individual 
X-ray s may be scattered multiple times. IMathis fc Lee I 
(1991) showed that for r sca > 1.3, there are more mul- 
tiply scattered photons than singly scattered. Mul- 
tiple scattering tends to broaden the halo in an 



energy-dependent fashion. The scattering cross sec- 
tion itself depends upon the X-ray en e rgy a nd the 
dust model; Table 1 of IMathis fc Leel l|1991|) shows 
that cr ;f , a . = 9.03 x lO" 23 ^" 2 , for d ust models 
such as IMathis Rumpl fc Nordsieckl l)1977t MRN77) or 
IDraine fc Led l|1984l) . For Afa» 4 x 10 22 cm~ 2 (see 
this corresponds to r sca = 3.6-E^y. Below 2 keV, we 
therefore expect both that the RG approximation will 
overestimate the total halo intensity and that the single 
scattering approximation will underestimate the radial 
extent of the halo. At 2 keV, however, the RG approxi- 
mation is adequate and although multiple scattering does 
not dominate, it will still observably broaden the halo. 
By 3 keV, however, single scattering entirely dominates. 

Explicitly calculating all possible scatterings is dif- 
ficult, but it is possible to calculate the radial pro- 
file including both si ngle and double scattering. B oth 
IMathis fc LeTI lITTiml) and IPredehl fc Klolel l(l99^ de- 
rived the general form for the total cross section including 
multiple dust scattering. Using the Gaussian approxima- 
tion to da/dU derived in ML91 in order to simplify the 
calculations, the singly-scattered profile is 
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where c\ = 1.1 cm 2 sr _1 , p is the grain density, in g cm~ 3 , 
a the grain size in pm, n(a) is the dust grain size dis- 
tribution, and E the X-ray energy in keV. This Gaus- 
sian approximation to the Rayleigh-Gans cross section is 
used only in this section to estimate the effect of multi- 
ple scattering; otherwise we use the full Rayleigh-Gans 
model with the form factor for spherical grains described 
by Equation [3 

Equation 0] can be extended using the recursion rela- 
tion given in Prcdehl & Klose ( 1996) to doubly scattered 
photons. This result holds for smoothly distributed dust 
between the xq and x\. In the following, we take ad- 
vantage of the small-angle nature of the scattering to 
simplify some of the trigonometry. As a result, the scat- 
tering over 9' does not include all possible values of 0' 
but is limited to 9\ « 20'. With that restriction, the 
profile for doubly-scattered photons is: 
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To our knowledge, this has no analytic solution. How- 
ever, it can be evaluated numerically if x% ^ 1 (otherwise 
the solution becomes numerically unstable, although the 
apparent singularity may only be due to the choice of 
coordinates). We considered smoothly distributed dust 



3 




2.5 keV: 
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Fig. 1. — Singly- and doubly-scattered surface brightness profiles through a column density of Nh = 4 x 10 22 cm -2 for 1.5 keV [Left] and 
2.5 keV [Right] X-rays, through a smoothly distributed MRN77 dust distribution from the observer to 90% of the distance to the source. 



between the observer (xo = 0) and 90% of the distance 
to the source (xi = 0.9). The scale of the effect can 
be seen in Figure ^ which shows the single and double- 
scattering profiles for Nh = 4 x 10 22 cm -2 at 1.5 and 2.5 
keV, assuming an MRN77 dust model. The total power 
in the doubly-scattered photons is not insignificant. The 
column density will produce an overall optical depth for 
scattering r = 1.25 at 1.5 keV and 0.42 at 2.5 keV. At 
1.5 keV, ~ 50% of the scattered p hotons will be multipl y 
scattered, and at 2.5 keV, ~ 20% (jMathis fc Lee 111991^ . 
However, the effect it has on the profile is significant only 
at large radii and lower energies. At 2.5 keV, even at 
large radii, the broadening effect is relatively weak. Due 
to the twin difficulties of calculating the double scatter- 
ing and the Mie cross sections we will focus on energies 
> 2.5 keV. 

3. DATA REDUCTION 

GX5-1 was observed by the Chandra ACIS-S for 7 ksec 
on August 6, 2000 (ObsID 109) and ~ 1.5 million counts 
were detected. Due to the design of the Chandra ACIS 
detector, even t "pileup" is a problem for bright sources 
lpavisll200l|) . For GX5-1, the pileup is so severe that 
no counts were detected within a 2.5" radius circle cen- 
tered at 18:01:08.217, -25:04:41.34 (J2000), after correct- 
ing the aspect solution following the Chandra website 
1 . We used a similar proced ure to that described in 
iSmith. Edgar fc Shaferl l|2002t SES02) to determine the 
source position for this observation, and estimate the er- 
ror in this procedure to be ~0.3" , in addition to the 0.6" 
(90% confidence) error in the Chandra aspect reconstruc- 
tion. Our result is 1.24" from the position determined 
for the IR counterpart (1 8:01:08.222, -25:04:42.58, J2000) 
by Uonker et all 1)2000)) and 0.7" from the ATCA po- 
sition (18:01:08.2 3 3, -25 :04:42.044, J2000) measured by 
iBerendsen et al. I ()2000|) . 

3.1. Pileup 

We now have to measure the radial extent of the pileup, 
to determine where we can begin to extract radial pro- 
files with confidence. One way to measure this effect is 
to examine the ratio of 'unphysical' grade 7 (which have 

1 http://cxc.harvard.edu/cal/ASPECT/fix_offset/fix_offset.cgi 




Fig. 2. — (a) Observed radial profile of the ratio of grade 7 events 
to grade 0, 2, 3, 4, and 6 events on ACIS-S3, in three bands: 1-2 
keV, 2-3 keV, and 3-4 keV. Near the source, pileup reduces the 
number of good events (defined as those with grades 0, 2, 3, 4, 
or 6) as multiple photons are treated as a single multi-pixel event, 
(b) Same, for grade to grades 2, 3, 4, and 6. For comparison 
purposes the average ratio measured in Cas A is shown as a flat 
line. 



detected charge in five or more pixels out of a 3x3 region) 
to the normally allowed grades 0,2,3,4 and 6 (which have 
charge in only one to four pixels) as a function of energy 
and angular distance from the source. As the pileup di- 
minishes, the ratio should approach a constant value if 
the counts are still dominated by the source, and then 
transition to the 'background' value far from the source, 
as shown in Figure EI a) for three energy bands. In Fig- 
ure[2Ib) we show the ratio of grade (single-pixel events) 
to grades 2,3,4, and 6 for the same three energy bands. 
In both cases, the ratios seem to indicate that within 
100" of the source, pileup is affecting the grade ratios. 
Beyond ~ 250" the ratio might change, possibly due to 
the increasing fraction of background photons. We lim- 
ited our analysis to radii beyond 100" to avoid possible 
contamination. 

3.2. Spectrum 

Measuring the spectrum of a strongly piled-up source 
is difficult, since the source counts cannot be directly ex- 
tracted. The so-called "transfer streak," created by mis- 
positioned counts that arrived during the frame transfer 
provide the only clean unpiled spectrum in the data. The 
transfer streak has not yet been explicitly calibrated by 
the Chandra X-ray Center team. Nonetheless an accu- 
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rate spectral measurement is crucial when subtracting 
the instrumental PSF as well as calculating the total 
dust column density. Our approach is similar to that 
described in SES02, except for the treatment of the back- 
ground. The source counts were extracted from a 8 pixel 
wide by 611 pixel long box that stretched across the 
CCD, with the exception of the region within 100" of 
GX5-1 (203 pixels), which was excluded to avoid pilcup. 
The allowed region contained 48465 counts, with a total 
observing time of 6868.7 seconds, or a total of 2146 3.2 
second CCD "frames." Each CCD frame, therefore, had 
~ 22.5 counts in the transfer streak, spread over a 8x611 
pixel box so pilcup in the readout streak is not signifi- 
cant. Since the 8 pixel wide box contains essentially all 
of the transfer streak, the total effective exposure time 
for the transfer streak is 2146 frames x 611 rows/frame 
x41 /is/row = 53.76 seconds. 

It is important to note that this is the exposure time for 
the source counts. The background flux (measured using 
boxes parallel but above and below the source boxes) is 
caused by scattered photons from GX5-1, cosmic rays, 
and other sources. These events could arrive during the 
transfer streak, but are much more likely to occur during 
the 3.2 s hold time, so the appropriate exposure time is 
the full 6868.7 seconds. This is an unusual situation: the 
exposure times for the source and background photons 
in the source spectrum are different, and none of the 
standard spectral tools (XSPEC or Sherpa) can explicitly 
deal with it. We dealt with this problem by treating both 
the source and background as having the full exposure 
time of 6868.7 seconds, and then scaling the source flux 
by the ratio of the exposure times, 6868.7/53.7495 = 
127.767. 

After extracting the source and background spec- 
tra, we found an adequate fit using an absorbed 
bremsstrahlung model (x 2 = 1.2) with TVjj = 4.0 x 
10 22 cm" 2 , kT = 10.5keV, and Fx (1-10 keV)= (5.2 ± 
0.1) x 10 -8 erg/cm 2 /s (absorption corrected). Fortu- 
itously, there was also a 9 ksec HETG observation of 
GX5-1 (Obsid 716) done July 18, 2000, only 19 days 
before the ACIS observation described here. Although 
GX5-1 is a variable source, the XTE All-Sky monitor 
shows little change between the two observations. In the 
1.3-3.0 keV channel, the RXTE count rate was 14.2±1.05 
during the HETG observation and 13. 4± 1.01 during our 
ACIS observation. In the 5.0-12.1 keV channel, the count 
rates were 32.3 ± 1.3 and 35.1 ± 1.3, respectively. 

The HETG observa tions of GX5-1 have been analyzed 
bv lUeda et al. 1 1)2005[) in order to measure the elemental 
absorption along the line of sight. As noted in §1, they 
found the total 7V H = 2.8±?;| x 10 22 cm~ 2 , based on ex- 
trapolation primarily from Mg, Si, S, and Fe absorption 
features. Their spectral model was more complicated 
than ours, but their total (absorption corrected) Fx(l- 
10 keV)= 4.3 x 10~ 8 erg/cm 2 /s. This is ~ 20% lower 
than our result from the transfer streak which may be 
due to true variation in the source. However, given the 
limited data available for the ACIS transfer streak, this 
20% variation may also be due to calibration uncertain- 
ties. In any event, although this overall flux uncertainty 
will be reflected in the measured dust column density, 
but does not affect the profile of the halo. 

3.3. Imaging Analysis 
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Fig. 3. — The CO emission observe d in the direction of GX5-1 
IDame. Hartmann fe Thaddeus 11200 II) . 



To calculate the radial profile of GX5-1, we first iden- 
tified the serendipitous sources from the image, using 
the CIAO routine celldetect. We then selected an en- 
ergy grid (1-4 keV, AE = 0.2 keV) and a radial grid 
(100 points from 10-800", log-spaced). We filtered the 
data using the measured CCD energy, and extracted 
the counts in the radial profile using concentric annuli 
centered on the source. After making exposure maps 
for each energy band, we filtered out the regions with 
serendipitous sources in both the data and each expo- 
sure map. We then extracted the effective "radial ex- 
posure" using the same concentric annuli, for each en- 
ergy band. Dividing the radial profile of the counts in 
each band by the effective area and exposure time gave 
the corrected radial profile in units of photons cm -2 s -1 
arcmin~ 2 per energy band. We then divided this by the 
observed source flux in each band to get the fractional 
radial profile in units of source fraction arcmin~ 2 . We fit 
this as the sum of the scattered halo, the Chandra PSF, 
and the background (both instrumental and cosmic). As 
with SES02, the Chandra PSF was measured from ob- 
servations of Her X-l. We also note that the Chandra 
PSF is small enough that we can neglect scattering of 
the X-ray halo by the telescope. Finally, although the 
cosmic background is vignetted while the instrumental 
background is not, we found that in these bandpasses we 
could simply fit them together as a flat continuum out 
to 800". 

3.4. CO and 21 cm observations 

The CO emission toward GX5-1 (Figure is con- 
fined to 4 broad lines that can be readily associated with 
three distinct regions of the Galaxy. First, the line peak- 
ing near -22 km s _1 almos t certainly arises from the 3 
kpc Expanding Arm (e.e. iBania Ill98fl() . one of the best- 
defined spiral arms in the Galaxy. It is seen tangentially 
at I ~ 24°, implying a Galactic radius of 3.46 kpc (as- 
suming a distance to the Galactic Center of 8.5 kpc). 
Since the inclination of the arm is unknown, we assumed 
it is circular to calculate its distance as 5.1 kpc at 5° lon- 
gitude. A second region is associated with the line at 22 
km s _1 , which arises from the so-called molecular ring. 
This is the region of high molecular cloud density roughly 
half-way between the Sun and the Galactic center. The 
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TABLE 1 

Molecular Concentrations toward GX5-1 from CO observations 



Cloud Velocity W co 2xN(H 2 ) N(Hl) Distance from Sun 

(km/s) (K km/s) (10 20 cm 2 ) (10 2Q cm 2 ) (kpc) 

Molecular Ring 22 40.4 146 34 4.7 (or 3.3; see JO) 

3 kpc Arm -22 24.4 88 (or 18*) 34 5.1 

Galactic Center (Tangent) region 184, 208 21.7 16+ 11 8.5 

Total Line of Sight 86.5 250 (or 180) 83 



+We ha ve reduced the expec ted H2 by a factor of 5, since molecular gas near the Galactic Center is thought to be anomalously bright in 
COfe.g. lSodroski et alll995h . 



near and far kinematic distances of this gas are 4.7 kpc 
and 12.1 kpc, with uncertainties of ±0.5 kpc imposed by 
the cloud-cloud velocity dispersion. At the far distance 
the molecular gas would lie 215 pc above the plane, more 
than 4 times the vertic al dispersion of molecular clouds 
({Bronfman et al."l ll988) , so the near distance is far more 
likely. 

It is worth noting that most of the molecular ring emis- 
sion closer to the plane at I ~ 5° is centered roughly 
10 km s _1 lower than that toward GX5-1; this (mildly) 
suggests that the molecular-ring emission toward GX5-1 
might be anomalously high by that amount. If so, its 
kinematic distance should be reduced to ~ 3.3 kpc; this 
possibility is also noted in Tabled We conclude that the 
molecular ring gas probably lies 3-5 kpc from the Sun. 
Given the extreme velocity crowding at this low longi- 
tude, it's possible that numerous clouds are spread over 
this range. A third region is associated with the other 
two partially-blended lines. These lines are at such high 
velocities (~ 200 km s _1 ) that the emitting region must 
lie near the tangent point, at a Galactic radius of 0.7 
kpc and thus near the Galactic center at a distance of 
~ 8.5 kpc. Distances for the main molecular concentra- 
tions along the line of sight to GX5-1 are summarized in 
Table □ 

The molecular column density, N(H2), corresponding 
to each of the CO lines can be calculated from their 
velocity-integrated intensities, Wco- F° r the molec- 
ular ring emission we use a standard Galactic value 
for X = NfHsl/Wm = 1-8 x 10 20 cm^K^km^s 
ijDame. Hartmann fc Thaddeus I l2001[) . Since there 
is evidence from both diffuse gamma ray emission 
I Blitz et al. I I1985D and far infrared dust omission 
I Sodroski et al. 1993) that molecular clouds in the Galac- 
tic center region are anomalously bright in CO, we re- 
duced A by a factor of 5 for the two high- velocity lines 
near the Galactic center. The molecular cloud in the 3 
kpc Arm might also be anomalously bright in CO, given 
the arm's proximity to the Galactic center and large 
non-circular motion. Our results are given in Table ^ 
along with the corresponding at omic column densities 
derive d from the 21 cm survey of H artmann fc Burton I 
(1997). The interpolated 21 cm spectrum toward the 
source shows velocity components similar to those seen 
in CO. We found a total H I column density along the 
lin e of sight of 0.83 x 10 22 cm~ 2 , in good agreement with 
the lDickev fc Lockmanl l|199(D value of 0.91 x 10 22 cm" 2 . 
The total gas column density, 2xN(H 2 ) + N(H I), is 
3.3 x 10 22 cm~ 2 , or 2.6 x 10 22 cm -2 , depending upon the 



value of X used for the 3 kpc Arm gas (see Table 

As with many X-ray binaries, the true distance 
to GX5-1 is not well determined. We use the 
iChristian fc Swank I l)1997|) upper limit of 9 kpc, based 
on its flux and calculated Eddington luminosity. Com- 
paring this upper limit with Table^shows that GX5-1 is 
almost certainly behind the 3 kpc arm and the molecular 
ring, but may be in front of, behind, or embedded in the 
Galactic center region. The radial profile of the X-ray 
halo depends on the dust position relative to the source 
flux, but since X-ray halos are due to forward scattering, 
dust behind the source is unimportant. 

4. RESULTS AND DISCUSSION 

GX5-1 is one of the brightest persistent X-ray sources 
in the Galaxy, and its relatively large absorption column 
density means that it has a substantial X-ray halo. The 
CO and H I observations also constrain the fit param- 
eters. As a result, GX5-1 is one of the best sources 
to use in testing IS dust models. We therefore col- 
lected a wide range of propos ed dust models: MRN77, 
iWeingartner fc Draine 1 120011 WD01), and the 15 mod- 
els proposed in lZubko. Dwek fc Arendt I l|2004 ZDA04) 
(most notably their BARE-GR-B model), as well as 
th e model including ext r emely large grains discussed 
in iWitt. Smith fc Dwek I l|20(m WSD01). These mod- 
els were developed to match optical, UV, and IR ob- 
servations of dust primarily found in diffuse clouds (or, 
in the case of WSD01, in-situ data) without considera- 
tion of their X-ray scattering properties. Figure HfLeft] 
shows the size distributions for three models, while Fig- 
ure ^jRight] shows the size distributions weighted by a 4 , 
making them proportional to the total X-ray scattering 
cross section following Equation|21 (a fter i n tegrat ing over 
all scattering angles; see iMathis fc Leel (jl99ljn . This 
shows that the differences in the largest grains must dom- 
inate the results. By systematically comparing all of 
these to the GX5-1 data, we can determine which models 
agree with the observed X-ray scattering, and which do 
not. However, it is important to realize that the dust 
grains scattering the X-rays from GX5-1 are primarily 
in dense clouds, and may have quite different properties 
from the dust in diffuse clouds. Therefore, this study 
cannot globally exclude dust models, but simply deter- 
mine which could describe the dust in dense clouds and 
which cannot. Our results are additionally limited by 
the extreme pile-up, which obscures the true profile when 
9 < 100". Dust very near the source (e.g. the clouds in 
the Galactic center region) will lead to scattered photons 
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in our insensitive "near region" (see also Figurc lTU)l In ad- 
dition, the number of scattered photons at every angle is 
proportional to FxNh, so errors in the measured flux lin- 
early correlate with errors in the dust column density. Fi- 
nally, the total column density tow ards the source is not 
well-k nown. The measurement by iChristian fe Swank I 
l|1997D of AT H = 2.54 ± 0.19 x 10 22 cm" 2 is reasonably 
precise, but its accuracy depends entirely on the spec- 
tral model assumed. The same is true of our best-fit 
value N K = 4.0 x I0 22 cm" 2 (see with the ad- 

ditional caveat that calibr ating a tran s fer str eak spec- 
trum is challenging. The lUeda et a l. ( 2005) result of 
7Vh= 2.8^1'g x 10 22 cm~ 2 , based on absorption features 
seen in high-resolution spectra, is accurate but far less 
precise. Considering the uncertainties on the Hi large- 
scale optical depth the CO-to-H2 conversion factor, the 
value of Nh = 3.3 x 10 22 cm -2 derived from the gas trac- 
ers is probably only good to within a factor of 2. Based 
on all these results, the most we can say with confidence 
is that ./Vh is likely in the range 2 — 4 x I0 22 cm -2 . If our 
X-ray halo results are outside this range, it would imply 
is that the dust/gas mass ratio along this line of sight 
is abnormally high or low. Therefore, we cannot impart 
too much significance to the overall normalization of our 
results. However, we can require that the best-fit column 
density be independent of energy, and that the shape of 
the profile conform to the observations. 

4.1. Smoothly Distributed Dust 

Although the CO data strongly implies the dust along 
the line of sight is clumped, we begin with the simplest X- 
ray halo model, where the dust is assumed to be smoothly 
distributed along the line of sight ( f(x) = 1), and that 
it sca tters the X-rays at most once ijPredehl fc Sch mitt 
1995). We begin by showing the radial profiles fit with 
some of the most common or best-fitting dust models, 
to show the quality of the results. Figure El shows the 
radial profile of GX5-1 at 1.5 ± 0.1 and 2.5 ± 0.1 keV, fit 
using a single-scattering smooth dust distribution. We 
show three different dust models: the classic MRN77 
model as well as the more recent WD01 (with Ry = 
3.1,6c = 6 x 10~ 5 ) and the BARE-GR-B model from 
ZDA04. The best-fit column densities were (respectively) 



2.9, 2.1, and 4.2 x 10 22 cm~ 2 at 1.5 keV and 2.7, 2.2, and 
3.5 x 10 22 cm- 2 at 2.5 keV. At 1.5 keV all the fits were 
inadequate, although the ZDA04 model was the best 
fit. Although at 1.5 keV Mie effects could affect the 
halo slightly, a fit including these did not significantly 
change the results. However, as noted in i)2.1l compar- 
ing with Figure n shows that including multiple scatter- 
ing might have improved the ZDA04 model fit. Unfortu- 
nately, including multiple scattering from smoothly dis- 
tributed dust in these fits is not yet implemented. How- 
ever, in the case of GX5-1 the dust is almost certainly 
not smoothly distributed and so this enhancement is left 
to a future paper. At 2.5 keV, Figure QlRight] showed 
that multiple scattering is far less important, and we 
see in Figure OjRight] that the ZDA04 model fit the ra- 
dial profile of GX5-1 reasonably well- x 2 .(ZDA04) = 1.8. 
The MRN77 and WD01 models were still poor fits with 
X 2 (MRN77) = 9.4 and x^(WD01) = 30. They both 
overestimated the scattering at small angles and under- 
estimated it at large angles far more than the ZDA04 
model. We then fit the GX5-1 radial profile data for 
energies between 2.5-3.5 keV (in 0.2 keV wide bins) to 
a model consisting of the Chandra PSF, a flat back- 
ground, and the scattered halo calculated for smoothly- 
distributed dust following each of the above dust models. 
In FigureEJwe show the xt values for each best-fit model. 
While none of the models are formally acceptable, clearly 
some are preferred. For example, none of the composite 
models of ZDA04 (those beginning with COMP-), with 
the possible exception of COMP-GR-B, are even close 
to a reduced xt ~ 1; while the BARE-GR-B model is 
a reasonably good fit {xt — 1-5), as would be expected 
from the results in Figure |3| 

While instructive, Figure is not conclusive. First, 
GX5-l's brightness means that the errors are probably 
not statistical but are instead dominated by calibration 
errors in the PSF and effective area. Second, the as- 
sumption of smoothly distributed dust disagrees strongly 
with the CO observations and is almost certainly in- 
correct. At the same time, the consistently poor fits 
found when using the composite models of ZDA04 sug- 
gests that these models poorly describe the observations. 
In Figure CfLeft] we show the fits to the radial profiles 
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Fig. 5. — [Left] Fractional halo profile of GX5-1 at 1.5 keV, fit with a (single-scattering) smooth distribution using the MRN77 (solid 
line), WD01 (dotted line), and ZDA04 (dashed line) dust models. Multiple scattering is not included; including it would increase the model 
scattering at large angles and decrease it at small angles. The power-law fit to the Chandra PSF is also shown as a solid line. [Right] Same, 
at 2.5 keV, except that multiple scattering is much less significant. 
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Dust Mode 

Fig. 6. — xi values for each dust model used, assuming smoothly- 
distributed dust along the line of sight and using the data between 
2.5-3.5 keV in 0.2 keV bins. At these energies, multiple scattering 
should be largely insignificant. 



using the COMP-GR-B, COMP-AC-S, and COMP-NC- 
FG models. All of these over-predict the scattering at 
small angles and under-predict it at large angles. Since 
larger grains create smaller halos, this suggests that these 
models have relatively too many large grains. Interest- 
ingly, however, the best-fit column densities (shown in 
Figure EfRight]) are almost all in the 2 — 4 x 10 22 cm -2 
range, with the exception of the last six models. In gen- 
eral, then, these models have problems not with the total 
mass of dust required, but rather the relative distribution 
of large and small grains. 

4.2. Clumpy Cloud Models 

Although the good smoothly-distributed ZDA04 
BARE-GR-B model fit could be taken as evidence for 
that model over the others considered, we must also con- 
sider the possibility that the match is at least partially 
fortuitous as the dust is almost certainly clumped into 
the same clouds seen in CO emission. In this case the 
halo profile depends upon the relative positions of the 
clouds and the source. We start by considering the sim- 



plest possible model, a single cloud with variable relative 
position and column density. The parameters were al- 
lowed to vary freely without any assumptions taken from 
the CO observations. This model again includes only sin- 
gle scattering. We begin with fits to the same three dust 
models (see Figure |SJ), and again find the best results 
from the ZDA04 BARE-GR-B model. At 2.5 keV, all 
models appears to fit equally well but only the ZDA04 
BARE-GR-B model has plausible fit parameters. The 
best-fit column densities for the MRN77, WD01, and 
BARE-GR-B models were 1.7,1.2,2.5 x 10 22 cm" 2 , re- 
spectively and the best-fit dust positions (in units of frac- 
tions of the total distance to the source) were 0.16, 0, and 
0.32. In both cases, the MRN77 and WD01 parameters 
are well outside the expected values. In addition, the 
BARE-GR-B model's best-fit column density appeared 
to be nearly constant with energy, unlike the other two 
models, and it had the lowest overall xt values. The best- 
fit ZDA04 BARE-GR-B model had N a = 2.5 x 10 22 cm" 2 
and position ~ 30% of the distance to GX5-1, with 
slight variation between 25 — 35% as a function of en- 
ergy. Nonetheless, even this fit is formally unacceptable 
(xt > 2) at most energies . Note that all fits strongly re - 
jected the result found bv lXiang. Zhang, fc Yaol (2005) , 
who found > 90% of the dust to be very near the source. 

We now consider a more complex spatial distribution of 
dust clouds, based on the CO data. The single cloud re- 
sults foreshadow possible consistency problems between 
the scattered X-rays and a "naive" interpretation of the 
CO data. As Tabled shows, multiple interpretations of 
the CO data itself are possible. The 3 kpc Arm may 
have substantial or insignificant amounts of dust, and 
the Molecular Ring dust may be at either 4.7 or 3.3 kpc. 
We considered all combinations of these options, and de- 
cided to use a value of 1.8 x 10 21 cm -2 for the 3 kpc Arm 
column density and a distance of 3.3 kpc for the Molec- 
ular Ring, as other choices led to a poorer fit to the data 
for all dust models. Of course, this may be an indication 
that there is a problem either in the calibration or the 
dust models, but we first examine if the closer distance 
can lead to a good result. 

In Figure [UJLeft] we show the halo profile fit using 
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Fig. 7.— [Left] Radial profile of GX5-1 at 2.5 kcV, with the best-fit ZDA04 COMP-GR-B, COMP-AC-S, and COMP-NC-FG models 
assuming smoothly-distributed dust. At large scattering angles these models all under-predict the total scattering, while at smaller angles 
(< f 50") the models overestimate the observed halo. [Right] The best-fit column density N^j for each model; the horizontal lines mark the 
expected upper and lower values. Most models find a reasonable total dust model, despite an overall poor fit. 
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Fig. 8. — Best-fit scattered X-ray halo to the GX5-f data for 
a single cloud at 2.5 keV. With the ZDA04 BARE-GR-B model, 
jVji= 2.5 X 10 22 and the cloud position is 33% of the distance to the 
source; \l = 2.1. The MRN77 and WD01 models, despite their 
similar appearance, are worse fits with \i = 3.7, 5.7 respectively. 



these preferred values from the CO observations to fix 
the cloud column densities and positions, and show the 
expected halos from a ZDA04 model considering single 
scattering only with three different positions. In this 
case the only variable fit parameter is the source flux. 
For these fits, the best-fit value is within 25% of the flux 
measured from the transfer streak, within the calibration 
uncertainty. Once again, the ZDA04 model fits were su- 
perior to either the MRN77 or WD01 models, and so 
for clarity only the ZDA04 model is shown. A distance 
of 15 kpc is clearly a poor fit, over-predicting the halo 
brightness at nearly all angles (a distance of ~ 15 kpc 
for GX5-1 would also imply the source has a consis- 
tently super-Eddington luminosity l)Christian fc Swank I 
1997)). The best-fit appears to be bracketed between the 
7-9 kpc distances, although neither of these is truly an 
acceptable fit. Only single-scattering was considered for 
these plots; multiple scattering would tend to broaden 
the halo slightly and might improve the fit slightly for 
the 7 kpc distance. However, as shown in Fig ureUlRight], 



the improvement would only be marginal and is unlikely 
to explain the entire difference. 

In Figure|5lRight] we again show the xt values for each 
best-fit model, this time for the multiple cloud model as 
suggested by the CO data. We assumed a distance for 
GX5-1 of 8 kpc, which gave the best overall results. As 
with Figure none of the models are formally accept- 
able, although some of the same models are again pre- 
ferred, most notably the ZDA04 BARE-GR-B model. As 
with the smoothly-distributed dust model, the compos- 
ite dust models of ZDA04, especially the "COMP-NC" 
models which use No Carbon, give very poor fits. 

5. CONCLUSIONS 

We have extracted the spectrum and scattered halo 
from GX5-1, one of the brightest Galactic X-ray bina- 
ries. Although the spectrum had to be extracted from 
the poorly-calibrated transfer streak, we were able to 
confirm the result due to a fortuitous HETG observation 
taken nearly contemporaneously with the ACIS observa- 
tion along with the RXTE All-sky Monitor data. GX5-1 
is so bright that pileup was evident out to a radius of 
100". Nonetheless, we were able to extract reliable sur- 
face brightness profiles between 100 — 800". GX5-1 has 
a large (but rather poorly determined) column density 
which is likely in the range 2 — 4 x 10 22 cm~ 2 . It is also 
close enough to the Galactic plane that velocity-resolved 
CO measurements exist which show the presence of a 
number of dense clouds along the line of sight. 

We have shown that GX5-l's column density is large 
enough that X-rays with energies below 2 keV are often 
scattered multiple times by dust along the line of sight. 
Since this is a complication to fitting, wc limited our fit- 
ting to energies above 2.5 keV where multiple scattering 
is much less significant. We compared our observations to 
a wide range of dust grain models, including for the first 
time the many models described in ZDA04. We found 
reasonable fits for some dust models assuming either a 
smooth dust distribution or a single "cloud" along the 
line of sight. No dust grain model fit the default param- 
eters determined from the CO data for any distance to 
GX5-1. However, the fit parameters found when using a 
single cloud suggested a re-examination of the CO data. 
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Fig. 9. — [Left] Observed radial profile (divided by source flux) at 2.5 keV with ZDA04 BARE-GR-B dust model using column densities 
and positions derived from the most plausible interpretation of the CO data with 3 possible source positions. [Right] \i values for fits to 
the data between 2.5-3.5 keV (in 0.2 keV bins) for each dust model, assuming the same parameters and a source distance of 8 kpc. 



The 3 kpc Arm (at 5.1 kpc) is on the boundary between 
the Galactic center, w here the CO is anomalously bright 
ijSodroski et al.lfl995^ . and the Galactic disk. If the gas 
in the 3 kpc Arm is in fact similar to the Galactic center 
emission, it would have only l/5th of the molecular gas 
(and associated dust) as predicted in Table H In this 
case, most of the dust would be the one in the Molecular 
Ring, which at 3.3 kpc is roughly 1/3 of the distance to 
the Galactic center and which has a predicted hydrogen 
column density of ~ 1.8 x 10 22 cm~ 2 . When using these 
values, we found better fits but somewhat surprisingly 
still not as good as the smooth dust distribution results. 
More data from X-ray binaries at or near the Galactic 
center are needed to see if in fact this interpretation of 
the CO data holds when compared to other X-ray halo 
observations in this area, or is merely special pleading. 

We also note our stron g disag reement with the results 
of iXiang. Zhang, fc Yao 1 1)2005|) regarding the dust den- 
sity along the line of sight. Although we are not as sen- 
sitive to dust near the source as their observation, the 
effects would still be noticeable in our results. We do 
not know the origin of this discrepancy, but note that 
they found significant amounts of dust near all of their 
sources. It seems possible that there may be calibration 
issues, especially in the zero-order of the HETG, which 
are in some part responsible for their result. 

We also found that some dust grain models strongly 
disagreed with the observations in all cases. Most no- 
tably, nearly all of the Composite dust models of ZDA04 
(models beginning with COMP-) disagreed strongly with 
the data (the sole exception is the COMP-GR-B model). 
The composite "No Carbon" (COMP-NC) models in par- 
ticular were uniformly poor descriptions. The BARE- 
GR-B model, which combines bare graphite and silicate 
grains with PAHs using B star abundances, was the over- 
all best fitting model in all cases. However, since no dust 
model or spatial distribution led to \t ~ 1 fits, and given 
remaining uncertainties in the calibration and the model- 
ing we can not state that this dust model is "preferred." 
The idea of directly measuring the IS dust mass, com- 
position, and size distribution is enticing enough that it 
appears in every proposal to observe X-ray halos. How- 
ever the history of astronomical X-ray scattering halo 




Fig. 10. — Predicted near-source X-ray halos at 2.5 keV using 
the ZDA04 BARE-GR-B model with total jV H = 3.91 x 10 22 cm" 2 . 
Three different dust position distributions are shown: smoothly 
distributed dust, a single cloud 1/3 of the distance to GX5-1, and 
three clouds at the positions predicted by the CO data, assuming 
a source distance of 8.51 kpc. 



results shows that there are difficulties arising from the 
many assumptions needed before the data can be mod- 
eled. Although not conclusive, our results show that 
adding CO and Hi observations do help in reducing the 
allowed model space and putting better constraints on 
the dust models. However, the uncertainties in the CO 
and Hi analysis must be included to sensibly use this ad- 
ditional data. We plan to test our results by analyzing 
more Galactic Center sources to see if, for example, they 
also imply that the 3 kpc Arm is CO-bright. Another 
way to break the existing near-degeneracy between dust 
models would be to obtain measurements of the radial 
profile between 10" — 100". Figure ITUI shows that varia- 
tions in dust grain positions affect the surface brightness 
at these angular distances far more than in the 100 — 800" 
region measurable with our data. We have just obtained 
Chandra HRC-I observations of GX5-1 which will mea- 
sure the radial profile in this range (albeit without energy 
resolution). 
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